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Abstract 

In a recent work, M.Kohout (M.Kohout, Int.J.Quant.Chem. 87, 12 2002) raised the important 
question of how to make a correct use of Bohm's approach for defining a quantum potential. In 
this work, by taking into account Kohout's results, we propose a general self-consistent iterative 
procedure for solving this problem. 
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I. INTRODUCTION 



Bohm 

mm 



im's formulation of quantum mechanics in terms of single-particle trajectories 



) has been, and still is, a continuous matter of dispute, often in rather philosophical 
terms above all concerning the fundamental meaning of 

quantum mechanics. At a concrete level such a theory has been used as a theoretical tool for 
understanding and interpreting several processes in different fields, from molecular ph ysic s to 



plasma physics, from scattering theory to simulation of quantum wires 
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2l|, to name a few. Recently M.Kohout j2^ raised the important question of how to treat the 
more rigorous and realistic 3A^-dimensional formulation of the quantum potential generated 
by electrons. Rigorously speaking, in an A^-particle system Bohm's potential is a 3A^- 
dimensional function and the one-particle Bohm's potential, usually considered in literature, 
is the simplest 3-dimensional reduced form which systematically does not take into account 
the effects due to the presence of other particles. He proposed a formal interpretation 
of the wavefunction as a product of a one-particle marginal function and a conditional 
many-particle function where somehow the effects due to the other particles are taken into 
account, then a formal expression of the quantum potential is obtained. However, as the 
author underlines, for the conditional many-particle function an explicit expression 
is required and this is a rather difficult problem . In this work, by taking into 
account Kohout results we build a reasonable initial guess for the many-particle potential 
and, following this choice, we develop an iterative self-consistent procedure for obtaining 
a general (numerical) expression of the potential. We restrict our analysis to a spinless 
system, in any case there exists the possibility to extend the procedure to wavefunctions 
which explicitly consider the spin variables. To conclude, we must underline that the 
intention of this work is simply to show that the widely used Bohm potential (used in the 
approximation of one particle) can actually be reasonably treated in its true form of many- 
particle and we propose a method to do so. It is not our intention to prove that such 
a procedure is preferable to other quantum approaches, such as Hartree-Fock or Density 
Functional Theory, for determining general many-body effects in electronic systems. 
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II. A BRIEF ACCOUNT OF BOHM'S THEORY 



Bohm's formulation of quantum mechanics in terms of single-particle trajectory is based 
on the assumption that the wavefunction determines the dynamics of more fundamental 
variables (hidden variables). The essence of the theory can be summarized by quoting 
Bohm's original work " The first step in developing this interpretation in a more explicit 
way is to associate with each electron a particle having precisely definable and continuously 
varying values of position and momentum" X]- Within this assumption, the system is 
described by its wavefunction '?/'(xi..., xjv, t) and the positions of its particles x^; these two 
fundamental quantities are governed respectively by Schrodinger equation 

= Hi; (1) 

where H, as usual, is the Hamiltonian of the system, and by the dynamical equation: 

d 

^x. = V. (2) 

where v,- = —Im ^ X^^ - This procedure leads to a non- Newtonian equation of motion 
which becomes Newtonian in the classical limit of ^ — * 0; the connection of such an ap- 
proach with quantum mechanics is provided by the fact that the quantum formalism auto- 
matically emerges from Bohm mechanics. In simple terms this approach " ...implies however 
the particle moves under the action of a force which is not entirely derivable from the clas- 
sical potential^ ^(x), but which also obtains contributions from the quantum mechanical 
potential..." 1]. Given the A^-particle electronic wavefunction whose general form can be 
written as tp(Il,t) = x(R, t)e*'^'^*^^'*^; R G 3?'^^, substituting it into the time-dependent 
Schrodinger equation and then separating real and imaginary part, once the velocity is 
defined as v = ~/ = ;;^V5', and p = IV^P = one obtains the equations: 

^p + V-(pv) = (3) 

and 

dS [V5(R,t)]2 



^ , 2^ -[mt) + g(R,t)]. (4) 

where m is the mass of the particle, 1^(R, t) is the potential characterizing the system 
(e.g. electrostatic for interacting fermions), (5(R, t) = —^ ^J^ is the Bohm potential and 
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V = Y2iLi gradients of the A^-particles. We restrict our analysis to the 

stationary case, thus we can write in terms of the wavefunction phase factor: 



2m 



[V{R) + Q{R)] . (5) 



III. KOHOUT'S FORMULATION OF THE MULTI-PARTICLE POTENTIAL 
PROBLEM 

Bohm's potential in three dimensions is used in many applications, within this framework 
p(r), a one-particle electron density, is defined as: 

p{r) = N |?/'(r,r2,r3 rN)prfr2rfr3....rfrN (6) 

{Q domain of definition of the system in real space) and the single-particle wavefunction 
writes: 

The wavefunction's form of Eq|3 corresponds to the procedure of separating Bohm's dynam- 
ical equations into independent and indistinguishable single-particle equations and describe 
a set of identical particles moving in an average potential where specific mutual interactions 
are neglected. In simple terms it is sufficient to describe only one particle, embedded in 
an average potential generated by the other particles, in order to automatically describe 

the whole system. In this case S'(ri, tn) = s(ri) + ....s^tm)- It must be noticed that, 

within the approximation done, what we called the one-particle Bohm potential, does 
not correspond to the Bohm potential for a system composed by only one par- 
ticle. In fact, if this was the case, Bohm equations could not be defined at the nodes of 
the electron wavefunction, instead the definition of ip{r) given in Eq|7|with p(r) defined by 
Eqini suggests that it is unlikely to have a ip which contains zeros at least for dense systems. 
The one-particle potential is the simplest approximation that can be used for describing an 
electronic system. Recently M.Kohout proposed another approach, by redefining the wave- 
function ?/;(r, r') (r' indicates the remaining — 1 particles of an N -particle system when 
one focuses the attention on particle r): 

^(r, r') = x(r, rY^^^''''''^ = 0(r) ■ P{r'\r)e"'^^''''''\ (8) 
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0(r) is the marginal function while /3(r'|r) is the conditional function, conditionally depen- 
dent on the coordinate r. These two functions satisfy the following normalization conditions: 



/ 

Jn 



(r)pdr=l (9) 

Jn 
and 

|/3(r'|r)|W = l;Vr. (10) 



/ 



IQ. 

By using this factorization, he obtains a formal expression of the potential Q(r, r') as a sum 
of the one-particle Bohm potential and a conditional multi-particle potential Qcond (r|r') = 
~ 2m^^r'|r) Yli=i N^i l^i^' \^) ' "^^^^ cxprcssiou barcs a rather difficult problem: to find a 
reasonable expression for /?(r'|r). In this work we will circumvent this problem by using 
a suitable expression for x(r, r') which recovers the properties of the factorized wavefunction 
introduced by Kohout, and leads to an iterative procedure for Q{r,r'); this is reported in 
the next sections. 



IV. MANY-PARTICLE DENSITY AND BOHM'S WAVEFUNCTION. 

As anticipated in previous sections our final aim is to develop a procedure for obtaining a 
Bohm's potential where many-particle effects are somehow incorporated. For this purpose, 
we redefine the wavefunction of the system by extending the form ■ip{r) — 0(r)e'^^'"^ for 
one-particle to ■i/)(ri, r2...rM) = ^(ri, r2...rM)e*'^^''^''"^ ' ""'^\ an M-particle wavefunction in 
an A^-particlc system, with M < N. We should also require to be antisymmetric and S 
to be symmetric with respect to the Ml possible pair permutations of the M particles; this 
requirement will preserve the antisymmetry of "ip. Later we will show that the symmetry 
of S corresponds to the fundamental physical property of indistinguishable particles. In 
the next sections, to remind the analogy with the one-particle case we will identify r with 
ri, thus 0(ri, r2...rM) = 0(r, r2...rM)- This basically means that given N particles, our 
system is characterized by (or alternatively, we are interested in considering) M-particle 
interactions which produce observable effects on the average behavior of the system, thus M 
particles must be treated explicitly; we can describe those effects by considering the M-th 
approximation, where M = 1, ...TV, the case M = 1 is the trivial one-particle case, M — 2 
counts two-particle effects etc.etc. For the moment we are not interested in ^(ri, r2...rM) 
whose role will be clear later on. Let us focus on (/)(r, r2...rM)- In the one-particle case 
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N ■ |0(r)p = p(r) = N J^(n~i) |^(r,r2,r3 rN)p(ir2(ir3....(irN = p(r), i.e. the average 

electron density of indistinguishable particles projected on the real (3- dimensional) space. 
In analogy we can define an M-particle electron density as : 

iV|0(r....rM)P = p(r....rM) = / |'0(r, r2..rM, Tm+i rN)prfrM+i....c?rN (11) 

where, in analogy to the one-particle case, one has: 

|0(r....rM)|'rfr....(irM = l. (12) 



This form of 0(r, r2...rM) satisfies the requirements of the factorization in marginal and 
conditional part as defined by Kohout: 

^{v, r2, ....rM) = 0(r) ■ '^(^'^^-^m) ^.,5(r,......M). (13) 

0(r) 

In fact the conditional function of EqlHl (3{r'\r), can be written as: /3(r'|r) = '^'''"'^(r)'"'^'' ^ind 
(as can be easily verified) one obtains J^^ |/9(r'|r)p(i*^r' = 1; Vr e Q, where r' = r2, r^, •••rM- 
This kind of factorization was already considered by Hunter 2^ who, once more, underlines 
that the nature of 0(r), as a marginal probability density function, makes unlikely the 
existence of zeros^ In terms of probability '^'''"'^(r)'''^'' is interpreted as the square root of 
the classical expression for the conditional probability density. The conditional probability 
density, i.e. the probability density for M — 1 particles given the position of the particle r, 

is written as the probability density of the particles r, r2 Fm, (|0(i", r2....rM)P), divided 

by the probability density of particle r, (|0(r)p). In general the mathematical structure 
of quantum mechanics leads to a non-commutative probability theory which coincides with 
the classical one only in case we treat commuting spaces !24]. Here for a commuting space 
we intend a set of variables representing physical quantities, such as positions, which do 
commute. In the language of quantum mechanics this means for example that position 
operators of the M-particles commute, i.e. a measurement of the position of particle 1 does 
not infiuence the measurement of the position of particle 2. The same example does not 
hold in the space of spins, however since we have restricted our analysis to spinless systems, 
we can apply the rules of classical probability. In R.ef . [2^ is reported a study about the 
possibility of defining a quantum correction which takes into account a non- commutative 
probability for the conditional function. Once we have formally defined the many-electron 
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wavefunction via Eq ll31 we can determine a reasonable initial guess for Bohm's potential, 
by defining a reasonable 0(r, r2....r]vi) : 

Q(r,r,....rM) = -^ J2 ^Um'^^- (14) 

V. BOHM'S DYNAMICAL SYSTEM FOR M NON-INDEPENDENT PARTICLES 

In the previous section we defined, in analogy to the one-particle case, an M-particle 
wavefunction for an A^-particle system. Clearly the larger the value of M the more difficult 
the determination of the wavefunction, although in many cases two or three-particle effects 
may be enough for a basic understanding of some physical properties. So far by defining 
the many-particle wavefunction we simply gave a first approximation for the potential Q, 
within, for example. Slater orbitals or similar appropriate approximations (e.g. plane-waves 
expansion ). The idea for an effective coupling of M electrons is to insert an 
approximation of Q, as an initial guess, in a system of equations which do couple 
the M electrons. This set of equations is obtained by inserting the M-particle wavefunction 
in Schrodinger equation and repeating the procedure followed for the one-particle case within 
Bohm's framework as it is shown next. In fact, by inserting ip{r,r2, •••rM) in Schrodinger 
equation for the stationary case, similarly to Eq|Sl one obtains: 

2^ 2^ = -y(r, ...rM) - Q(r, ...rM) (15) 

where — ^(r, ...rM) is the electrostatic potential experienced by each single particle. In this 
form, EqElis not very useful thus we should find a way for a simplification which can be 
based on physical well founded hypothesis. The first thing to notice is that Vr.-S'(r, ...rM) = 
Vr.(r, ...rM) represents, within Bohm's framework, the velocity field of particle r^ which 
depends upon (or alternatively, is influenced by) the positions of the other M — 1 particles 
as well. We can also notice that — ^(r, ...rM) can be represented as the electrostatic potential 
per particle and can be expressed in good approximation as: 

nr,...rM)= E ^ 

i=l,M 

In EqUni p(w) is the one-particle electron density, obtained by integrating the M-particle 
density over M — 1 variables. The meaning of Eq^Jis that in a system of indistinguishable 
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p(w) 



-dw 



w 



(16) 



particles the average electrostatic potential experienced by one-particle is the same as that 
experienced by another particle. The "average" character is obtained by considering the 
one-particle electron density instead of the M-particle one. The use of a simpler rather than 
a a more complicated expression for V can be justified by the fact that we know how to 
express electrostatic properties and we know that Eq^J is reasonable. What is unknown 
are the quantum effects represented by Q, for this reason we can take known quantities in 
their simplest approximation as far as they are known to be reasonable. In any case the 
important fact is that following Eq]16lone can write: 

V{r, ..fm) = V{r) + Vir^) + .... + V{rM)- (17) 

The expression for ^(r, ..fm) given in EqElis particularly useful for simplifying EqllSl In 
fact EqJ15lcan be written as : 

(18) 

where ]g(5(r, r2..., tm) is the quantum potential per interacting particle. At this point we 
can proceed having in mind the following points: 

EqJlHIis a 3M-dimensional non linear partial differential equation and we would like to have a 
solution which somehow recovers the dynamical physics contained in Bohm's approach. The 
natural simplification of Eq^J is to decompose it in a system of M one-particle equations 
coupled through ^(r, r2...rM) and Q{r, r2..., rm), where each equation describes the squared 
modulus of the velocity field of a particular electron Vi6'(r, ...fm) subject to a potential 
per particle [V{ri) + ]g(5(r, r2..., tm)] • Clearly the sum of solutions of single equations is 
also a solution for the initial one. In mathematical terms this kind of solution represents 
only one possible way (a particular solution) , while in physical terms this represents the 
direct extension of Bohm's one-particle dynamical description to a system of M-particle non- 
independent electrons . In this way we obtain a set of equations coupling the M interacting 
electrons as we postulated at the beginning of this section. The M-equation system takes the 
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|V.g(r..,.rM)|- 

2m 



^ — 1 A/T 



V{ri) + — Q(r,r2...,rM) 



form: 

(1) [v.(r,...rM)]^ = 1^i%^m)1! = _ [v{r) + -^g(r, r,..., Tm)] 

(2) K(r,...rM)]^ = '"^""t"^"^'' = - [Vir,) + jjQir,r,...,rM)] 

(3) 

(4) 



(M)K,(r,...rM)]^ = '^"-"t-""^' = - [^(^m) + jjQ{r,r,...,r^)] . 

In the next section we show that, by using EqUHl a self-consistent procedure for Q{r, ...fm) 
can be determined. 

VI. SELF-CONSISTENT ITERATIVE PROCEDURE FOR Q(r, . ..fm) 

The basic idea is to start from the approximative form of Q = Qo(r, •••rM) obtained 
via the M-particle density within a Slater orbitals or similar approach and then intro- 
duce it in equation (1) of system ^\ Next by solving the first order non linear par- 
tial differential equation for 5'(r, ...fm) with respect to r one find an "initial" value for 
S{r, ...Fm) = 5'i(r, ...fm)- At this point 5'i(r, ...tm) is substituted in equation (2) and leads 
to a new value of Q, Qi{r, ...tm) = ^^''^"^^^m "'"'^''^ ~^ [^(^2)]- The next step consists in substi- 
tuting Qi in equation (3) and solve it (as for equation (1) ) with respect to to find a new 
S'(r, ...tm) = S2{r, ...tm), the procedure is repeated until the convergence of Q meets a given 
criterion of acceptance. It must be noticed that the procedure is iterative in the sense that 
one should go through all the M-equations of the system and then take the final Q or S for 
repeating the process by going through the M equations again as is illustrated for the case 

M = 3 and M = 4 in Fig^ Finally Q{r, ri fm) can be integrated with respect to M — 1 

variable and so be reduced to a one particle potential where the many-particles effects are 
integrated out (averaged). However, in practical terms this procedure is not straightforward 
since a first order non linear partial differential equations must be solved, for example with 
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respect to r = (x, y, z) at the very first step: 



1 
2m 



9a; J \ dy J \ dz 



V{r) + ^Qo{r,...rM) 



(20) 

and formally a solution of the following form should be obtained: 

S{r, ...fm) = /(r, r2, ...rivi) + G{r2, ...t-m) + const. (21) 

G{r2, ••.rM) is a function constant with respect to r and const is a constant which can be 
neglected since we are interested in gradients of S. G(r2, ...rM) cannot be obtained by solving 
a single equation since in principle one should solve the whole system of equations and find 
the global solution (basically G'(r2, ...fm) would be an artifact of the procedure). On the 
other hand to find a global solution would represent a formidable, if not impossible, task. 
We need a well founded approximation which allow us to remove the artifact represented by 
the function G(r2, ...tm). At this point we notice that there exists a fundamental physical 

property of the system which can reasonably solve this problem. In fact S{r, ri tm) must 

be symmetric under any permutation of the particles. This fundamental property is a 
direct consequence of the fact that the particles are indistinguishable and implies that the 
velocity field of one particle can be obtained from that of another via particle permutation 
(see Appendix). 

The symmetry of the total wavefunction can be preserved by taking (p{r,r-j_, ....Fm) anti- 
symmetric with respect to any pair exchange of the M particles. The property of symmetry 
for S allows us to make the iterative procedure possible and physically reasonable. In fact 
we solve Eql2()l with respect to r and obtain S{r) as: 

^i(r,...rM) =/(r,r2,...rM) (22) 

formally neglecting, for the moment, the part G{r2, ...rM) of Eql21l Then in order to obtain 
a global solution which satisfies the permutation (symmetry) criterion we write the solution 
for S at the first iterative step as : 

Sitot{r, r2, ...rM) = ^ PtSi{r, ...Tm) (23) 

i 

where Yli Pi is the sum over all the possible M-coordinate pair permutations and Pi is the 
permutation operator for the i-th permutation. This procedure, by making a fully symmetric 

10 



solution S, does not require to consider the term G{r2, •••I'm) which would formally come 
from a direct mathematical approach. This point can be considered from an alternative point 
of view ; we extract a solution of a multi- variable function problem by considering only one 
equation instead of considering all the equations. Such a solution is valid for the variable 
r (integration variable) as long as the other variables are considered as constants (fixed at 
a particular parametric value). Next, by using the physical condition of indistinguishable 
particles, we can extend, reasonably well, this solution to the whole set of variable by making 
the global solution a fully symmetric function of the whole set of coordinates. We can 
then go to the next equation to find Q, make it fully symmetric with the same procedure 
we used before for S and go to the next equation, again solve for S, make it symmetric 
and proceed on in the same way. A formal argument for the symmetry of S is given in 
Appendix. Technically the major problem is the fact that in general to solve a non linear 
partial differential equation in three dimensions is not an easy task, however Eq|2ni is the 
well known equation of the eikonal in an anisotropic media. In this case the boundary 
condition for S must be assigned when the Schrodinger equation is formulated and should 
depend on the particular nature of the system under consideration. The problem of the 
solution of this equation goes beyond the purposes of this work, however there exists a 



massive list of references dealing with this prob 



numerical techniques (see for example 
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em for both mathematical properties and 
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0,13,0 



and references 



therein). In particular the fast sweeping algorithms of R,ef. j.Sflj l . and the robust algorithms 
for multidimensional Hamilton- Jacobi equations of Refs.l3il.l32.l33l represent an extremely 
useful approach to a computer implementation of this procedure. In particular these latter 
algorithms, although dealing with the more general non stationary problem, in principle can 
be adapted to the stationary case and in general would make it possible also to extend the 
procedure to non stationary cases. 



VII. CONCLUSIONS 



We propose a general method to treat Bohm's potential in its 3A^-dimensional rigorous 
form. Inevitably, there are several physical approximations that one should accept; we use a 
spinless system, thus we can apply the concept of classical probability to a quantum system 
in a stationary state. This choice offers technical advantages but is not physically obvi- 
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ous; in case the spin variables are explicitly considered, the basis of our procedure remains 
valid, but we should find an opportune way to define a more general wavefunction where the 
conditional function is determined by non-commutative probability principles as suggested 



in Ref. 



2J]. The separation of the many-electron equation in single but mutually depen- 



dent one-particle equations is not unique from a rigorous mathematical point of view but 
it is based on a reasonable physical approximation. From a technical-mathematical point 
of view the major problem is the solution of the non-linear partial differential equations 
for S (eikonal equation). As it is commented in several textbooks of quantum mechanics 
(see for example 0]), to solve this equation is not an easy task, however the literature 
relative to the solution of the eikonal equation is sufficiently large and several methods 
are available. For small M (e.g. 2,3) it should be possible to apply the procedure with 
a reasonable effort; moreover there exists systems which can be described with simplified 
unidimensional models such as electrons in a one dimensional wire or an isolated atom 
considered spherically symmetric. Other simple examples are those where it is possible to 
proceed to a separation of variables (i.e./(xi, Xj, yi, yj, Zi, zj) = X{xi, Xj)+Y{yi, yj) + Z{zi, zj) 
or f{xi, Xj,yi, yj, zi, Zj) = X{xi, Xj)Y{yi, yj)Z{zi, zj)) in this case the eikonal equation is re- 
duced to independent ordinary first order differential equations which can be easily solved. 
The simple cases listed above can be used as a first approximation of more complicated sys- 
tems. As stated before, it is important to underline that the intention of this work is that of 
proposing a procedure to properly treat and use Bohm potential in fields where it is currently 
(and extensively) employed for practical applications. We do not claim that this procedure is 
computationally or methodologically more convenient than others in solving general many- 
electron problems , however, this method may also represent a complementary theoretical 
approach to the standard ones used in current research; due to the deterministic interpreta- 
tion of Quantum Mechanics, on which the method is based, effects which cannot be describe 
by standard Hartree-Fock or DFT may be revealed. This issue, anyway, involves a much 
deeper analysis which goes beyond the purpose of this work and will be possibly treated 
elsewhere. 
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IX. APPENDIX 

For simplicity let us consider a two-particle system (r, r'). Let us suppose that from Eg QUI 
for (r,r'), integrated with respect to r, the following solution is found: 

/(r,r') = (?(r,r') + Mr)+/(r) (24) 

where g{r,r') = g{r',r) and l{r') is unknown since it represents the constant (with respect 
to r) obtained by solving the equation with respect to r. From EqEU it follows that the 
formal expression for the velocity field for the particle r is 

Vr(r,r') = Vr[(?(r,r') + /i(r)] (25) 

this means 

Vr(r,rO =Vi(r,rO+V2(r) (26) 

while for particle r' we have: 

Vr,(r,r') = Vr'[^7(r,r') + /(r)] = Vi(r', r) + V3(r'). (27) 

Because the particles are indistinguishable one must obtain Vr'(r', r) from Vr(r, r') by permu- 
tation of (r, r') thus the unknown function Z(r') must coincide with h{r'). As one can easily 
see, by taking /(r, r') = g{r, r') + h{r) and making it symmetric with respect to the permu- 
tation (r, r') (r', r), we reach the same result for l{r'). In the most general case g{r, r') is 
not necessarily equal to g{r', r), however the same procedure (and principles) can be applied 
with the only difference that g must be also made symmetric, i.e. gfin{r,r') = ''"^ 
and hfin{r, r') = ^iilitMLl^ case of more than 2 particles the same argument can be used 
by taking into account all possible particle permutations. 
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FIG. 1: The figure represents the pictorial iUustration of the self-consistent iterative procedure 
shown in this work. For STEP we intent the process of going through all the 3 (4) equations. Next 
by using the final Q or S, obtained at the first STEP, it is possible to start a second STEP. Sij 
(Qij) is the function S calculated at different stages of the procedure, i corresponds to the number 

of equations solved within the system , i = corresponds to the starting (initial guess) expression 
for Q, j indicates the number of global iterations (STEP). 
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